Geachte heer Biesheuvel, 
beste Simon, 


Dank voor je leuke artikel in Euclides 94-3 en de uitdaging die je aan het eind van het artikel 
het land instuurt. Als (mede-)auteur van het telhoofdstuk in Moderne Wiskunde (1led) kan 
ik het natuurlijk niet laten om op je uitnodiging in te gaan. 

Over opgave 6 wil ik het niet hebben. Die is gewoon fout uitgewerkt. Maar aan opgave 
43 hangt een ander verhaal. Het formuleren van zo’n opgave is niet zo moeilijk. En als je 
niet direct ook een uitwerking maakt, kun je je enorm vergissen in de moeilijkheidsgraad. En 
dat is hier gebeurd. Helaas heeft de uitwerkingenmaker de opgave niet goed geïnterpreteerd 
en daardoor niet aan de bel getrokken, met als gevolg dat er een veel te lastige opdracht in 
het boek terecht is gekomen. 

Tot zover de geschiedenis. Maar daar kan het natuurlijk niet bij blijven. Je schrijft dat 
je een antwoord hebt gevonden en dat je (redelijk) zeker bent van je resultaat. Mijn hoogste 
doel is om die latente onzekerheid weg te nemen. Ik probeer op twee manieren het juiste 
antwoord te produceren. De tweede manier laat zich ook redelijk eenvoudig generaliseren. 


1 Formulering van het telprobleem 


De opgave in het boek is geformuleerd in termen van een geordend rijtje kleuren. Maar 
voor het vervolg is het makkelijker om te praten over getallen. We vragen naar het aantal 


verschillende permutaties van de getallen 1, 1, 2, 2, .…., n, n waarin twee opeenvolgende 
elementen verschillend zijn. Zo’n permutatie heet (ben ik deze vakantie achter gekomen) een 
Carlitz-permutatie van de multiset (1, 1, 2, 2, .…., n, n). Laten we dit aantal maar noteren 


met cp(n). De opgave vraagt dus om cp(6) te bepalen. 


2 Methode 1: Brute force 


Als nadenken te lastig of te foutgevoelig is, dan wil het schrijven van een programmaatje nog 
wel eens uitkomst bieden. Met behulp van een redelijk standaard backtracking algoritme zijn 
alle Carlitz-permutaties te genereren. Tellen is dan eenvoudig geworden. Het resultaat staat 
in de tabel van figuur 1. 
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Figuur 1: Computerresultaten 


Dit zou een einde moeten maken aan je eventuele onzekerheid. Of het slaat een deukje in 
je zelfvertrouwen, maar dat kan ik me nauwelijks voorstellen. 


3 Encyclopedie 


Ken je The On-Line Encyclopedia of Integer Sequences (https://oeis.org)? Een geweldige 
site. Als je het beginstuk van een rij hebt, dan levert deze site je vast een heleboel informatie 
over de rest van de rij, over contexten waarin deze rij opduikt, directe formules en recurrente 
betrekkingen. Het feit dat deze site de getallen uit figuur 1 herkent en weet voort te zetten 
geeft mij de overtuiging dat mijn uitwerking van het algoritme correct is. 


4 Methode 2: redeneren 


Een programma levert het juiste antwoord. Maar je wordt er niet veel slimmer van. Sowieso 
vraagt het genereren van alle (Carlitz-)permutaties voor wat grotere waarden van n veel te 
veel rekentijd. n = 6 is nog wel te doen, maar veel groter moet het (op mijn machine) niet 
worden. 

Het valt (mij althans) niet mee om in één keer het aantal Carlitz permutaties te berekenen. 
Maar het lukt wel met een aantal tussenstappen. 

Het is niet moeilijk om het aantal (gewone) permutaties van de multiset 
(1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6) te berekenen. Dat gaat min of meer standaard: 


12! 


2e 


Maar daar zitten ook permutaties bij waar opeenvolgende elementen wel gelijk zijn. Ik noem 
twee gelijke elementen die naast elkaar staat een dubbel. Fen Carlitz-permutatie is dus een 
permutatie met 0 dubbels. We beginnen met het makkelijkste geval (6 dubbels) en werken 
toe naar O dubbels. Als we daar zijn, dan zijn we er. 

Laat nu p(k) het aantal permutaties van (1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6) zijn waarbij 
de eerste k cijfers dubbels vormen en de andere cijfers lossen. Bijvoorbeeld p(4) is het aantal 
Carlitz-permutaties van de acht symbolen (11, 22, 33, 44, 5, 5, 6, 6). 

Merk op dat cp(6) = p(0). 


4.1 Zes dubbels 


Een permutatie van (1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6) met zes dubbels is in feite een 
(gewone) permutatie van de zes (verschillende) symbolen 11, 22, 33, 44, 55, 66. En dat zijn 
er 6! = 720. Dus 

p(6) = 720 (1) 


4.2 Vijf dubbels 


Met 6 als los getal berekenen we het aantal Carlitz-permutaties. We hebben nu de zeven 
| 


symbolen 11, 22, 33, 44, 55, 6, 6. Het aantal permutaties van deze symbolen is 5 =— 2520. 

Maar daar zitten ook de permutaties bij waar de twee zessen naast elkaar staan, ofwel de 

permutaties met zes dubbels. Maar dat aantal hadden we al berekend: p(6) = 6! = 720. Dus 
7! 


p(5) = zj — P(6) = 1800 @) 


4.3 Vier dubbels 


Met 5 en 6 als losse getallen berekenen we het aantal Carlitz-permutaties. We hebben nu 
de acht symbolen 11, 22, 33, 44, 5, 5, 6, 6. Het aantal permutaties van deze symbolen is 
eid = 10800. Maar daar zitten ook de permutaties bij waar de twee 6-en en/of de twee 5-en 


naast elkaar staan. Maar die aantallen hebben we al berekend: 


e twee 5-en naast elkaar; de 6-en los: p(5) = 1800 (zie 2). 

e twee 6-en naast elkaar; de 5-en los: dit zijn er natuurlijk net zoveel als de 5-en dubbel 
en de 6-en los, dus ook p(5) = 1800. 

e zowel de 5-en als de 6-en dubbel: p(6) = 720 (zie 1). 


Zo komen we op 8 
p(4) = PI — 2: p(5) — p(6) = 5760 (3) 


4.4 Drie dubbels 


Met 4, 5 en 6 als losse getallen berekenen we het aantal Carlitz-permutaties. We hebben nu 
de negen symbolen 11, 22, 33, 4, 4, 5, 5, 6, 6. Het aantal permutaties van deze symbolen 
is Ce Maar daar moeten weer een aantal vanaf: 


e de permutaties waarbij één van de cijfers 4, 5 of 6 een dubbel vormt (en de anderen 
niet). Dat zijn er (5) -p(4). 
e de permutaties waarbij twee van de cijfers 4, 5 of 6 een dubbel vormen (en de andere 


niet). Dat zijn er (5) p(5) 
e de permutaties waarbij de cijfers 4, 5 en 6 alle drie een dubbel vormen. Dat zijn er 


(3) - p(6) 


Zo komen we tot 


PB) == (1) (5) - (5) 0) — 21900 4) 


4.5 Twee dubbels 


Het patroon is inmiddels wel zo’n beetje duidelijk. Bij de tien symbolen 11, 22, 3, 3, 4, 4, 5, 5, 6, 6 
eerst maar weer alle permutaties: . En daar vanaf halen: 


e één van de vier symbolen 3, 4, 5 of 6 dubbel: () -p(3). 
e twee van de vier symbolen 3, 4, 5 of 6 dubbel: 6) -p(4). 
e drie van de vier symbolen 3, 4, 5 of 6 dubbel: 6) p(5). 
e alle vier de symbolen 3, 4, 5 en 6 dubbel: (5) -p(6). 


Samen 


4.6 Eén dubbel 


Maar in één keer 


(1) 20+ (5) + (5) vr (5) 0) + (Ì) vO) — ros 


4.7 Nul dubbels 


Dan nu waar het om gaat 


emt) = 0 = TE) (5) ve) vr (vr) (6-0 — 201000 
@) 


5 Generalisatie 


Het is aardig om op een systematische manier een opgave op te lossen, maar het wordt pas 
echt leuk als daaruit ook een algemeen patroon valt af te leiden, zodat ook grotere versies van 
het probleem effectief (en efficiënt) zijn door te rekenen. 

Laat nu 


p(n,k) 


het aantal Carlitz-permutaties van de multiset 
(11, 22, …, kk, k1, k+1,...,n, n) 


Dan voldoet dit aan de recurrente betrekkingen 


p(n‚n) =n! 
e voor 0 Sk <n: 
n—k 
(2n —k)! n—k . 
mere ? i p(n,‚k+ú) 


en 
cp(n) = p(n,0) 
Achteraf was het misschien handiger geweest om niet te kijken naar het aantal dubbelen, maar 


naar het aantal lossen. Dus laat q(n,t) = p(n‚n —t). Voor q hebben we dan de recurrente 
betrekkingen 


q(n,0) = nl! 


e voor 0 <t <n: 


q(n,t) 


p(n, n= t) 
= _{ vervang k door n —t } 


vl D lr) 
(2n — (n—t))! DS (TED). p(n, (n —t) +1) 
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cp(n) = q(n,n) 


Dat ziet er al aanmerkelijk werkbaarder uit. Dit is ook eenvoudig te implementeren. In 
pseudo-code: 


FUNCTION cp (int n) { 
ARRAY a [O .. nl; // alil = q[n,i 


al0] := n!; 
FOR i = 1 .. n DO { 
alij := (n + i)! / 2"i; 


FOR j =O .… ú - 1) DO { 
alil := alil — binom(i,j)*aljl; 
} 
ds 
RETURN aln]; 
} 


Een paar resultaten voor q(n,t) met n verticaal en t horizontaal: 

































































0 1 2 3 4 5) 6 7 8 
1 1 0 
2 2 1 2 
5) 6 6 12 30 
4 24 36 84 246 864 
5 120 240 660 2220 8760 39480 
6 720 1800 5760 21960 96480 478080 2631600 
7 5040 | 15120 | 55440 | 236880 | 1149120 | 6219360 \ 37079280 | 241133760 
8 | 40320 | 141120 | 584640 | 2772000 \ 14736960 | 86597280 | 556426080 | 3876561360 | 29083420800 





6 Alternatieven 


In het op de site van de Encyclopedia aangegeven artikel van Eriksson en Martin worden 
twee andere aanpakken beschreven. De inclusie-exclusie methode lijkt veel op wat hiervoor 
is beschreven. Ze doen het iets anders en komen dan tot een (gesloten) uitdrukking met 
een (ingewikkelde) sommatie. Uit eindelijk verschilt de uitwerking dan niet zoveel van mijn 
implementatie. 

De aanpak met geordende permutatie-patronen is uitermate elegant. Dat bedenk je (ik 
althans) niet zelf. 





